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1 Introduction 

Gravitational theories based on quadratic curvature (QC) invariants were first consid- 
ered long ago in the context related to gauge invariance in [fl] as well as in (2'| soon 
after the advent of general theory of relativity fl3] . Since then their study have been 
motivated by various aspects of general theory of relativity and cosmology. QC terms 
are also motivated by the effective action for the low energy limit of string theories 
E|. Although the QC terms lead to violation of unitarity, and ghost degree of free- 
dom, they are favorable in the renormalization in the perturbative quantum gravity 
at one loop level [5|. The gravitational models with QC terms also have desirable 
cosmological features such as they provide a framework for inflation [6| and with 
nonlinear curvature models it is possible to avoid cosmological singularity Q. Non- 
linear gravitational lagrangians also have the peculiar property that they can be cast 
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into Einstein theory interacting with additional matter fields [8 |. Moreover, in con- 
trast to the cosmological models which introduces dark energy and dark matter, var- 
ious nonlinear gravitational models provide alternatives to account for the observed 
late time acceleration of the Universe ||9l . 

Currently, there has been a further interest in the modified gravitational theories 
that based on the actions involving QC terms especially in three dimensions. They 
are motivated by the search for a consistent model of quantum theory of gravity. In 
particular, in three dimensions New massive gravity (NMG) ifTTl theory involves a 
particular combination of QC terms complementing the usual Einstein-Hilbert ac- 
tion. In contrast to the topologically massive gravity (TMG) iflOl which introduces 
mass by addition of a Chern-Simons term to Einstein-Hilbert action, at the linearized 
level, NMG lagrangian is equivalent to Pauli-Fierz free, massive, spin-2 lagrangian. 
In addition, NMG theory also has parity-preserving field equations and it leads to a 
ghost-free quantum theory which is unitary at the tree level. More recently, critical 
gravity models, that can be considered extension of NMG theory to four dimensions, 
lead to renormalizable gravity theory by inclusion of QC terms and a cosmological 
constant lfl2l . 

The NMG action has peculiar properties, such as the trace of the metric field equa- 
tions yields back the QC lagrangian density and therefore the trace leads to second 
order partial derivatives in the metric components relative to a natural coframe. This 
particular property is crucial in showing that the linearized equations are equivalent 
to the Pauli-Fierz massive lagrangian and hence introduction of mass to graviton. Re- 
cently, there are research on extension of this property to QC lagrangians in higher di- 
mensions B13U14L The gravitational lagrangians having field equations of sixth order 
partial derivatives of the metric components are investigated in order to single out la- 
grangian densities that are composed of cubic curvature terms as well as terms having 
square of the derivatives of the curvature tensor llT4l . In JT3 J 14] it has been reported 
that in four dimensions, it is possible to construct QC actions whose metric field 
equations are of second order and the trace of the metric field equation yielding back 
the quadratic part of the original action one started with. QC lagrangians having this 
particular property has been constructed explicitly by using three curvature-squared 
terms, namely scalar curvature squared, Ricci-squared and Gauss-Bonnet term. The 
resulting QC lagrangians contains two free-parameters. Though, the generalization 
of NMG lagrangian to higher dimensions will be presented in arbitrary dimensions 
n > 3 and the lagrangian densities whose equations of motion lead to second order 
trace are constructed using a different approach allowing spacetime dimensions n = 3, 
n = 4 and n > 4 to be studied separately with the emphasis on the common structure 
of the fourth order lagrange multiplier terms. In particular, in three dimensions, the 
orthonormal coframe formulation of field equations allows one to provide an inter- 
esting mathematical connection between the field equations of the TMG and NMG 
theories. 

The gravitational lagrangians considered in what follows are based on QC in- 
variants, in fact, represents a narrow sub-class of modified models which include 
functions of these QC invariants. In particular, for the modified QC lagrangian of the 
form f(R, G) where G is the Gauss-Bonnet scalar (defined below as G * 1 = in 
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n dimensions) see the reviews, for example [15] and [16|. In the following, it will 
be shown that, in general, such a modified QC model can be put into a form which 
is dynamically equivalent to original QC lagrangian coupled with additional scalar 
field(s). 

The outline of the paper is as follows. In the next section, the first order con- 
strained formalism is highlighted using the exterior forms relative to orthonormal 
coframe. This section also serves to fix the notation that will be used throughout. The 
exterior differential form notation and the geometrical conventions used are taken 
from H17II18I . see also the appendix in [20] for a concise introduction. The general 
formula presented are specialized to QC gravity models in the following two sections. 
In the subsequent sections, the gravitational lagrangians in dimensions n > 4 and 
n < 4 having second order trace are constructed with emphasis on structure of their 
Lagrange multiplier term which contribute to the metric (or equivalently coframe) 
equations, n > 4, n = 4 and n = 3 cases are discussed in the respective separate sec- 
tions. In the section devoted to three dimensions, a formal mathematical relation be- 
tween field equations of NMG and TMG are constructed in terms of the geometrical 
quantities. In particular, the conditions under which TMG is exact "square root" of 
the NMG is presented. 

2 Variational derivatives relative to an orthonormal coframe 

Before delving into the explicit mathematical expression, for convenience, basic grav- 
itational variables are defined briefly here. The basis coframe 1 -forms are denoted 
by 0" and the metric tensor takes the form g = rj a p9 a <g) 9@. The exterior product of 
the basis 1-forms will be abbreviated as 6 a A 0$ A ■ ■ ■ = 9 a P ". The covariant exterior 
derivative D acts on tensor-valued forms. The torsion 2-form a = ^T^^d^ can be 
defined as 

& a =D9 a = d9 cc + (O cc j3 A9 13 (1) 

where T a ^ are the components of the torsion tensor and co a p are connection 1-forms. 
Relative to an orthonormal coframe metric compatibility for the connection 1-forms 
reads co a p + (Qp a — 0. In terms of Curvature 2-forms the Cartan's second structure 
equations reads 

flp = \ R V V = da> a p + C0« A (0% (2) 

where R a p iJLV are the components of the Riemann tensor. The contractions of the 
forms and tensor-valued forms are defined with the help of the contraction operator 
i ea = i a where e a is basis frame fields that are metric dual to the basis 1-forms: 

iccQ^ = 8a- * is Hodge dual that defines an inner product for two p-forms in terms 
of the metric tensor. In terms of the Hodge dual, the oriented volume element reads 
*1 = Tj\g\dx° A • • • Aflfx'" -1 ' in n dimensions. Ricci 1-form can be defined as R a = 

ipQP a = R^p^d^ 1 whereas the scalar curvature R can be written as R = i a R a - All 
the other tensorial objects will be introduced as they required in terms of the basic 
quantities defined here. 
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In this section the scheme of calculation of variational derivative relative to an 
orthonormal coframe will briefly presented [20.21]. In later sections, the general for- 
mulae of this section will be applied to all considerations in what follows. As will be 
apparent in the following sections, the formulation of the QC metric field equations 
relative to an orthonormal coframe leads to compact form that has an advantage in 
the algebraic manipulations involved compared to the similar calculations relative to 
a coordinate basis, cf. for example, [22 1. 

The gravitational lagrangian densities that will be studied below are all of the 
form 

jgf(») = jS?(»)[0« a)«] (3) 

are «-forms in n dimensional Lorentzian manifold. The first order constrained formal- 
ism where coframe and connection 1 -forms are taken as independent gravitational 
variables 1 33 34- 1 will be used. All the lagrangians considered will therefore depend 
on the basic gravitational variables, coframe and connection 1-forms. The lagrangians 
will depend on exterior derivative of these variables, namely, d9 a and dco a p through 
the torsion and curvature forms respectively. 

The corresponding field equations for the Riemannian metric is deduced from 
the coframe variation of the lagrangian (|3), subject to the condition that the torsion 
tensor & a = D9 a vanishes. Relative to an orthonormal coframe where the metric 
components are constants, the metric compatibility of the connection 1-forms reads 
(O a p + (Op a =0. This algebraic condition can simply be imposed by antisymmetriz- 
ing the coefficients of 8(Q a p whereas the torsion-free condition, © a = which is a 
dynamical constraint on the gravitational variable 9 a , can be implemented into the 
resultant field equations by extending the original lagrangian density to include a 
Lagrange multipliers (n — 2)-forms X a as 

^' [n) [6 a ,CO a pM=^ {n) +^ a ^®a- (4) 

The total variational derivative of the extended lagrangian density with respect to 
independent variables has, up to a disregarded boundary term, the general form 

8&' {n) [e a ,(o a p,x a ] = 8co al3 A[n^-^(e a A^-e^AX a )] 

+8B a A(*& a +DX a ) + SX a A0 a (5) 

where the convenient definitions 

n ap^^_ *Sf« = ^— (6) 

8(Q a p 89 a 

have been introduced. It is also convenient to define (n — 1)— form*!?™ = SJrf'^ /89 a 
= ¥S a +DX a . It is symmetric as a result of the local Lorentz invariance of the ba- 
sis coframe: 9^ A *E a — 9 a A = and it is covariantly constant as a result of 
diffeomorphism invariance of lagrangian composed of curvature scalars: D * E a = 0. 

In the orthonormal coframe formulation, the metric field equations derived from 
coframe variational derivative, namely the vacuum equations *E a = are defined by 
©. The vacuum field equations *E a = involve the exterior covariant derivative of 
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the lagrange multiplier term which can be expressed in terms of the other independent 
fields by making use of the variational derivatives of the extended lagrangian with 
respect to the connection 1-forms. The lagrange multiplier term, X a is a vector-valued 
(n — 2)-form while TI a ^ is antisymmetric bivector-valued (n — 1 )-form and they have 
equal number of components ll34l . The equivalence can be expressed precisely as 

X? =2i a n a v + l -eV M^n^ (7) 

which can be obtained by calculating two successive contractions of equations de- 
rived from variational derivative of the extended lagrangian with respect to the con- 
nection subject to the constraint vanishing torsion constraint. & a = results from the 
variational derivative of the extended lagrangian with respect to the lagrange multi- 
plier (n — 2)-forms. Because the contribution of the Lagrange multiplier term to the 
metric field equations induced by the coframe variations is given by DX a , the proper- 
ties of the tensor- valued (n — 1) form n a P becomes crucial for the study of lagrangian 
densities below since it determines the fourth order terms. Hence, it is convenient to 
derive some of the general properties that will be used later. 

The contribution of the Lagrange multiplier to the trace of metric field equations 
is 6 a ADX a which can also be expressible in terms of the bivector valued form n a $ 
using the general expression for it as 

e a ADX a = -2Di p (e a Ani 3a ) (8) 

In particular, this important result implies that in the case in which n^ a satisfy 9p A 
TIP® = 0, the La grange multiplier terms do not contribute to the trace of the metric 
field equations. 

As another property, the invariance under the general coordinate transformations 
require D*E a = which involves D 2 X@ . This term can conveniently be expressed 
in terms of bivector-valued form TI a ^ as 

D 2 X p =2Q li a Ai^ a (9) 

where the first Bianchi identity, D 2 0^ = Q^ a A a = 0, has been used. 

The covariant exterior derivative of bivector-valued 1-form can be calculated us- 
ing the field equations for the connection 1-form. In a metric theory, D corresponds 
to covariant exterior derivative of torsion-free metric compatible connection 1-form 
and the covariant exterior derivative of the equations for the connection 1 -form yields 

Dn a P = ~(e a ADi p -e p adx") (io) 

after D9 a — @ a — has been used. It is interesting to note that the right hand side 
of ( [Tol l vanishes identically for the QC lagrangians having second order trace stud- 
ied below. The general scheme of calculations presented here is specialized to QC 
lagrangians in the following sections. 
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3 Field equations for QC models and their trace 

As will be apparent in what follows, the field equations for QC lagrangians depend on 
the dimensionality of the spacetime n and therefore the dimensions n > 4 and n = 4 
and n = 3 are to be studied separately. The presentation of the field equations for the 
QC lagrangians allows one to study the field equations in its most general form. It is 
convenient to start with n > 4 dimensions and then specialize to lower dimensions. 
In n > 4 spacetime dimensions, the most general QC gravitational lagrangian will be 
parameterized by the following individual QC lagrangians 

JSfW = a^ n) +b^ n) + c^ n) 

= aR 2 *l+bR a A*R a + cQ a pA*Q at3 (11) 

involving the linear combinations all three of the QC scalars with three coupling 
parameters a,b,c. (The numerals and the Latin subscripts for lagrangian densities are 
not tensorial and are used as discriminating labels.) Using the field equations given 
in the previous section, the equations for the QC lagrangian (fTTT i can be written in the 
form 

Q a p A P *X? P - fjSf W +DA/ 1 = (12) 
where DA/ 1 = dXj 1 + G)^ v A A r v and the bivector- valued 2-form X t aP is given by 

X t ap =2cQ al3 + e a A(bR 13 -aRO^-OP A(bR a -aRO a ). (13) 

Since the field equations (fT2l do not contain the number of dimensions n explicitly, 
they are, in fact, valid for all n. For the calculations in the following sections, the 
auxiliary anti-symmetric tensor-valued 2-form X a P is very convenient and it is also 
to be used in the calculation of lagrange multiplier term in terms of the gravitational 
field variables as well. For QC lagrangians, it has the general expression as the lin- 
ear combinations of curvature 2-forms and Ricci 1 -forms and the scalar curvature. 
Explicitly, it can equivalently be introduced by the defining relation 

n a P =D*X a l } = d*X a P + a a ^A*X^ + a^A*X a L l . (14) 

The general QC metric field equations (fl~2) then can be cast into the following re- 
markably concise form 

DX{ l + *Tf[R,R a ,n a P}=0 (15) 

which separates naturally the terms that are linear in curvature (the lagrange multi- 

ocB 

plier terms) and QC terms. In terms of the contraction of the auxiliary 2-form X t 
and with the curvature 2-form, (n — l)-form *T, can be defined as 

* T? [R,R a , £2°*] = ~(^n ap )A*X t aP + l -Q a ^Ai^xf (16) 
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for the total QC lagrangian in (fTTT i. The auxiliary (n — 1 )-form *T t ^ has a mathemati- 
cal structure that is formally analogous to energy-momentum form of a generic mini- 
mally coupled matter 2-form field, F — jF a p9 a P, expressed relative to an orthonor- 
mal coframe (see, e.g., H17II18I ) and henceforth will be called energy-momentum 
(n — 1 ) -form for the QC terms. 

In the general case corresponding to (1TTV the trace of the energy-momentum form 
^[R,^,^^] can be found as 

T, =T t ^[R 1 R a ,Q a P]*l = i(«-4)i2°^A*X, a/3 . (17) 

The general trace expression vanishes identically for n = 4 which in fact reflects the 
dimensional dependent nature of the general QC lagrangian (fTTI ). Moreover, note here 

(n) 

that S£f can conveniently be rewritten in terms of contractions of auxiliary 2-form 

aft 

X t H with curvature 2-form as 

i? f W = i*2 aj3 A*x/^ (18) 

also that the lagrange multiplier term in (fl5t which is linear in curvature and it com- 
prises all the terms that are fourth order in the partial derivatives of the metric com- 
ponents relative to a coordinate basis. Therefore, the form of the field equation ( fT3b 
suggest that the trace of the field equations is either fourth or second order in partial 
derivatives of the metric and the only way to have a second order trace is to have DA/ 1 
term to have vanishing trace. In the latter case, the trace becomes proportional to the 
general lagrangian density (fTTI as will be elaborated below. Finally, note that the co- 
variant exterior derivative (or equivalently the divergence) of the energy-momentum 
form can be calculated, by making use ofn" 13 =D*X t aP in the general formulae (0 
of the previous section. 

The addition of Einstein-Hilbert and/or cosmological terms is of considerable 
interest in applications and will be considered in especially in the sections devoted to 
three dimensions and the expression of the energy-momentum form for QC terms in 
terms of the auxiliary 2-form X a $ then allows one to write down the corresponding 
terms in the field equations. Explicitly, if an Einstein-Hilbert lagrangian density of 
the form jR * 1 is added to the general QC lagrangians (fTTI . then the corresponding 
field equations can effectively be obtained by the formal replacement 

x «P^ x «P + ^ e ap (19) 

(<7 and K are constants) in (fT2l . It this convenient to recall at this point that the Ein- 
stein (n— 1 )-form *G a — *{R a — jR9 a ) can be expressed in terms of the contraction 
*G a = —j£2 a p A P * 9 a P. Addition of n-form A * 1 to ( fTTI ). with A being a cosmo- 
logical constant, leads to the replacement 

*Tf ^*Tf+X*9^ (20) 

for the energy-momentum terms of the QC terms in the general field equations ( fT31 ). 
These formal replacements apply to any QC lagrangian in n > 3 dimensions studied 
below. 
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It is well known that the particular QC n— form, namely the Gauss-Bonnet term 

^ = \Q a pAQ liV A*d a ^ v (21) 

leads to second order field equations in the partial derivatives of the metric and 
consequently the corresponding trace of the metric equations is also second order 
fl9l . In four dimensions, Gauss-Bonnet the variational derivative vanishes identically 
whereas in three dimensions, the term itself vanishes identically. Although the Gauss- 
Bonnet QC term will not be considered below such a term can easily be included into 
the general analysis in n > 4 dimensions, for which the variational derivative with 
respect to connection 1-form leads to 

nf B =D{£l nv A*e a ^ v ). (22) 

and this expression vanishes identically as a result of second Bianchi identities. In 
turn, the corresponding lagrange multiplier (n — l)-form X£ B vanishes identically. 
Therefore the contribution of J^gb to coframe variational derivative with respect to 
80x is second order in the partial derivatives of the metric components and is of the 
form i*2 MV A Q afi A *9 a Pv^ . 

Next, it is convenient to the study the trace of QC field equations in the most 
general case, by wedging the metric field Eqn. (fT5t from the left by basis coframe 
1-form and using dTSb , it is possible to show that the trace of the metric equations 
simplifies to the expression 

- (n - 4-)Q a p A *X? P + M A DA/ 1 =0. (23) 

The Lagrange multiplier term can further be simplified by making use of the general 
formulae dS) with fl"^ specialized to that of derived from X"^ in (fT3l as 

OfiAXf 1 — 2i a D * {[(n — 2)b + 2c]R a + [2{n - l)a + b]R9 a } (24) 

where fl"' 3 +Tlf a = has been used. Consequently, the trace of the general field 
equations for QC lagrangian takes the form 

(n - 4)_Sf ( W - 2Di a D *{[(«- 2)b + 2c]R a + [2(n - 1 )a + b]R9 a } = 0. (25) 

As a result for all the QC gravitational actions, the number of spacetime dimensions 
n = 4 turns out to be critical dimension. In particular, for n = 4 the trace equations is 
simply the homogeneous wave equation for the scalar curvature in the most general 
case. The contribution of the Lagrange multiplier form involves the covariant exterior 
derivatives of terms linear in Ricci forms and scalar curvature. As mentioned above, 
this implies that there are in general two different cases regarding the trace in which 
it involves fourth and second order partial derivatives of the metric. 

In the fourth order case where the Lagrange multiplier term contributes to the 
trace, the trace can further be simplified by making use of the contracted second 
Bianchi identity D * G a = 0. In doing so, the trace then takes the form 

{n - 4)_Sf/ B) - [4(« - 1 )a + nb + 2c]d * dR = (26) 
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assuming that the coefficient of the d*dR term does not vanish identically. Even in 
this case, the number of spacetime dimension n = 4 is singled out with regard to 
the trace of the field equations for the most general QC lagrangian. For n = 4, (l26l i 
reduces tod*dR = 0, unless the coefficient in the square brackets is zero and becomes 
fourth order in the derivative of the metric components. The subcases for which the 
term in the square brackets vanishes identically will briefly be studied from slightly 
different point of view below. 

In the case where the trace is second order in the derivatives of the metric, the 
contribution of the Lagrange multiplier term can be made to vanish by tuning the 
coupling parameters a,b,c. Note also that from the result in ((26) that in this case the 
trace of the field equations yield back the lagrangian from which the field equations 
are derived up to a numerical constant depending on the dimension. In n > 4 space- 
time dimensions, there are three distinct subcases each of which corresponds to QC 
lagrangians that are inequivalent. These cases can be derived as linear combinations 
of the QC lagrangians as follows. 



(1) c = 0, -S?i2 = + bJzfj for which the metric field equations read 

Q aj3 A!" * [9 a A (bRP +aR9P)- 9? A{bR a + aR9 a )]-i^^ +£>Af 2 = (27) 



Consequently, using this auxiliary 2-form given in Q, it is possible to show the 
expression for Lagrange multiplier form can be simplified to 



Using the general considerations of the previous section, the trace of the metric 
equations can be found as 



(« - 4)^ - 2Di a D *{(n- 2)bR a + [2(« - 1 )a + b]R9 a } = 0. (30) 



In contrast to the fourth order trace-case above, instead of using the contracted 
second Bianchi identity to simplify the trace expression, in order for the trace to 
be of second order, one requires that the lagrange multiplier term to be propor- 
tional to the contracted Bianchi identity. This can only happen when the expres- 
sion inside the curly brackets above becomes proportional to Einstein 1-form. In 
other words, with the judicious choice of the constants 



the trace expression reduces to (n — 4) Jz?^ = as a result of D * G' 1 =0. This 
subcase also turns out to be also persist in three and four dimensions as well as 
will be presented below and therefore the result generalizes in fact to dimensions 

in) 

n > 3. It is worth noting that the QC lagrangian density _z 12 , w i tn tne coupling 



with the corresponding auxiliary bivector-valued forms given by 



n" 2 P = D* [9 a A (bRp+aROP) - 9 P A (bR a +aR9 a ) 



(28) 




(29) 
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constants satisfying OTl ). d29i l yields l a =< *DL a and that it is possible to rewrite 

in) 

QC lagrangian ^f 12 in the following equivalent forms 



rfiW 
-^12 



4(n-l 



-,6 = 1 



= -L a AL p A*9 afi = -G a A*L a (32) 



where L a =R a — ~, ,M 9 a is the Schouten 1-form, defined in n > 3 dimensions 

2(n-l) — 

11241 and G" =R a — |/?0 a is the Einstein 1-form. Various properties of (l32l such 
as the variational derivatives with respect to L a and G a have been studied in 
11251 in more general context of Riemann-Cartan geometry using the language of 
exterior differential forms. 

Finally, note that, with the addition of an appropriate cosmological constant to 
d32l . the lagrangian density becomes the critical theory [12| which has recently 
been introduced for n = 4. 

(2) b = 0, _Sf/" 3 = JSf/^ = a&[ n) + c^ 3 (n) has the metric equations 

2Q ap Ai^ * {cQ^+aRO^-P^f +£>Af 3 = (33) 
with the corresponding auxiliary bivector-valued forms given by 

n"f = 2D * (cQ a P+aRd afi ). (34) 
and thus the trace is now given by 

(n-4)«5f 1 ( 3 ) -2Di a D*[cR a + (n-\)aR9 a ] = 0. (35) 
Analogous to the first case, if the two non-zero constants are related by 



a 



c (n-1) 



(36) 



in) 

then trace expression again reduces to _£fj 3 — by means of the contracted sec- 
ond Bianchi identity D * G^ = 0. 
(3) a = 0, Jz?/ n) = JSf^ = 2>jSf 2 W + c_Sf 3 W . The metric field equations that follow 
from this density takes the form 

Q a p Af * [2cQ a P +b(0 a ARP - 9? AR a )]-i^^ ] +DX^ = (37) 

with the corresponding auxiliary bivector-valued forms given by 

n 2 °f =D*[2c£2 a P + b(9 a AR 13 -9? AR a )]. (38) 

As in the previous two cases, for the trace of the field equations to be of second 
order, the constants b,c are to be chosen such that the Lagrange multiplier terms 
add up to contracted second Bianchi identity. Explicitly, the trace of the metric 
field equations are 

(n - 4)Jz? 2 ( 3 " ] - 2Di a D * { [2c + (n - 2)b]R a + bR9 a } = 

Thus, if the constants b and c satisfy 2c + nb = 0, the fourth order Lagrange 
multiplier term vanishes identically. 
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4 Modified QC models 

In order to study modified QC models which involve some functions of QC terms 
relative to an orthonormal coframe, one first defines the relevant QC scalars Q,P,K 
by making use of ( fTTT i as 

S^ n) =Q*l, ^ n) =P*l, ^ n) =K*l. (39) 

Now, it is convenient to consider the gravitational lagrangian that depends on the 
quadratic curvature invariants Q,P,K in the form 

^w.=/G3,W*i (40) 

with the assumption that / satisfies some regularity condition stated below. Note that 
although it is possible to start with a function / depends also on other QC terms, 
such as, Gauss-Bonnet term and/or Weyl curvature-squared term. For the sake of the 
simplicity of the argument, such terms will not be considered. 

in) 

It is easy to write the total variational derivative of the lagrangian ScIJj in the 
form 

8 ^!"li. =fQ8^ n) +fp8^ n) +fK8^ n) + [f-Qf Q -Pf P -Kf K ]8*l. (41) 

The second term on the right hand side contributes to coframe equations only and in 
particular it is, by definition, the Legendre transform of the function f(Q,P,K). This 
identification then allows one to introduce the scalar fields 

Consequently, the variational derivative of the modified lagrangian JS^j in fiTT i be- 
comes equivalent to that of the Scalar- Tensor(ST) type lagrangian 

if P = frjSf/" 5 + 2 if 2 (n) + 3 ^ 3 (,,) + V(0i , h, <h) * 1 (43) 

provided that the potential term V(0i,02,03) f° r the scalar fields is taken to be the 
Legendre transform of f(Q,P,K) up to an irrelevant sign, namely 

V(h,h,k)=f-Qf Q -Pfp-Kf K . (44) 

Similar to the well-known f(R) theory case, a dependence of the function / on a 
particular QC invariant introduces a corresponding scalar field via an appropriate 
Legendre transform. The change of variables {/q , fp , /k } i-> { 01 , 02 1 03 } requires that 
the Legendre transform to be non-singular. The fulfilment of the requirement however 
depends on the explicit form of the function /. 

In the above general framework, the dynamically equivalent ST form of the field 
equations for the modified gravitational lagrangian (|4D1 can concisely be written as 

01 * Ef + </> 2 * + 3 * £ 3 a + V (01 , 02, 03) * e a 

+2d<j>i Ai p D*X" P +2D[(ipd<j) l ) *X" P ] 
+2d<h AipD*X" P +2D[(i p d<h) *xf /3 ] 
+2dfoAi(tD*X" P +2D[(i p d<j)3)*X" 13 } =0 (45) 
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where *E" = are the field equations corresponding to the QC lagrangian densi- 
ties Jz?/"^ for k = 1,2,3 respectively and x"^ are the corresponding bi- vector valued 
forms defined in the previous section. 

The trace of the metric field equations (l45l l in this general case now determines 
the equation that scalar fields satisfy. However, the dynamical nature of the scalar 
fields depends crucially on the explicit form of the function f(Q,P,K). 

Finally, in the light of the above discussion, it is worth to point out that mod- 
ification of the QC topological term -Sfgg = G * 1 to the form GB = /(G) * 1 
inevitably induces a scalar field and consequently, it is possible to write a dynamically 
equivalent lagrangian density in the form 0-2gg + V (<p ) * 1 in terms of the scalar field 
(j) = fc defined as in the above examples. The elaboration of ideas related to the ST 
equivalents of such modified QC models will be taken up elsewhere. 

5 n > 4 dimensions 

Although the form of the field equations derived above are the same for all dimen- 
sions, apparently, the number of spacetime dimensions turns out to be a crucial pa- 
rameter for the field equations and consequently in the calculation of their trace. Note 
that in n > 4 dimensions the QC lagrangian densities jSf^, %n and Jzf^ are not 
equivalent since they can not be related by an addition of a total derivative QC term. 
The results of the cases (l)-(3) allow one to construct more general QC actions. 

An arbitrary linear combination of the densities -S?/^' ', and -S^" w ' tn three 

arbitrary constant constitutes the most general QC lagrangian density with second 
order trace, or to put it more precisely, the trace is proportional the lagrangian density 
from which the field equations are derived (compare the construction, for example, 
given in HP ). As a result of the above analysis of QC lagrangians, it is easy to 
construct the most general QC lagrangian having second order trace in terms of J2?i2, 
jSfi3 and Jz?23 with three free parameters. In the most general case, the total lagrangian 
density takes the form of the linear combination 



with p,q,r being the free parameters of the pure QC lagrangian. Note that the field 
equations for JCg are given by 



,(«) 




g 



nr)Q a p A *Q 



(46) 



Q aj3 A ^ * Xf ~ i»^ n) + DXg = 



(47) 



where 



nr)Q aP + 6 a A (p + r)RP 




fl^A (p + r)R 



(48) 
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and Tig = D * X g . Note that the trace of the metric equations imply that when the 
metric field equations are satisfied the action vanishes identically since the tracing the 

(n) 

metric equations gives J2g ' = for n 9= 4. This result can in fact be used to further 
simplify the metric equations (|47| |. 

In three and four dimensions, on the other hand, only two of the cases out of 
total (l)-(3) are independent and one of the cases can be written in terms of the 
remaining two. Note also that, in all the cases, the existence of a second order trace, 
or equivalently, the vanishing of the lagrange multiplier terms, requires this term to 
involve the Einstein form, in contrast, in three dimensions it involves Schouten form. 
On the other hand, neither the case (2) nor (3) seems to have been investigated so far 
as the QC lagrangian that generalizes the second order trace property of the NMG 
lagrangian to dimensions n > 4. 

As another important example of QC lagrangians in n > 4 dimensions, which 
will also be of interest in four dimension is the Weyl tensor- squared QC lagrangian 
density. From the point of view adopted above, it corresponds to the particular values 
of the parameters p,q,r. Weyl 2-form C a $ — jC a ^ ^ v 0^ v can be defined in terms of 
curvature 2-forms and its contractions by ifPTl 

c afi = Qa p _ _!_^a AR p _ Q p AR a } + J_ (49) 

(n — 2) \n—\)(n — 2) 

for n > 4 dimensions. By contracting this equation side by side one finds that 

~4 n) = C ap A *C a P = Q ap A *i2 - j^Ra A *R a + -, ^ -J? * 1 (50) 

which is suitable for the study of the field equations in the above framework. The 
corresponding metric field equations can be written down immediately as 

DX& + *T£=0 (51) 

with —2* C a P . In this case, the lagrange multiplier term turns out to have a 
simple expression in terms of Cotton 2-form as a result of the identity 

n a P = 2(n-3) (0 « A ^_ « A ^ ) (52) 
(n-2) 

for the bivector-valued form TI W . Consequently, the use of variational derivative 
with respect to connection 1-form in the general formulae (0 immediately yields the 
vector-valued lagrange multiplier (n — 2)-form as 



K = T^*C a (53) 
[n — 2 



where C a = DL a Cotton 2-form and L a =R a - ~, 1 ,, R0 a is the Schouten 1-form in 

2{n— 1) 

n dimensions introduced above. The expression ((53) can be found, by first obtaining 
DC"! 3 = ^(0 a AC p -9 p AC") (54) 
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by taking the covariant exterior derivative of (l49l and then using the duality relation 
*C a P = \C^ V A *e a £ v together with the fact that D * 9 a ^ v = (26). Explicitly, in 
terms of (|52] >. It is worth stressing at this point that C a = DL a in the expression (1531 
is four dimensional Cotton 2-form and have remarkable resemblance to its three di- 
mensional analogue, namely to the corresponding expression for the NMG theory 
cf. Eqn. (l64l below. Using the properties of the Weyl 2-form, the energy-momentum 
forms *7m for Weyl 2-form can also be put into the following convenient and sim- 
plified form. Consequently, using all these result, it can be shown that (l50l l can be 
rewritten in the form 

A Jf^D*C»--^—R a A*C» a + *Tg=Q (55) 
(n-2) (n-2) 

where is defined as 

* T£ ee - {f C afi ) A *C afi + C ap A * C ap . (56) 

Finally, note that for the lagrangian rjx-^w + R*1, the field equations can be 
written in the form 

( *D *D - ^%m 2 \ R a = — — r * {R<x A *C» a ) + (-iy+"T£ (57) 
\ 2{n — 3) ) n — 3 

by taking the Hodge dual of the equation. The sign of the last term on the right hand 
side in d57l i results from the identity ** = (— \) s+p ( n ~P>id with square of the Hodge 
dual acting on a p— form in n dimensional space with the metric signature s. 



6 Four dimensions 



In four dimensions n = 4, owing to the fact that the variational derivative of the Gauss- 
Bonnet lagrangian vanishes identically, 8J?Qg = 0, the variational derivatives of the 
terms R 2 * 1, R a A *R a and Q a p A *Q a P are not all independent [27 1 . As a result, the 
most general form of the QC lagrangians in four dimensions comprises simply one 
of the linear combinations (i,k = 1,2,3) considered above. 

It is convenient to consider the following relation among the QC lagrangians hav- 
ing second order trace 



«.(»)_ «(») _ _2(n-3) 
GB w ~ (n-2) 



R" A *Rrr — ■ 



R z *l 



4(n-l) 

which holds for n > 4. In four dimensions, this implies that 

^ = R a A *R a - l -R 2 * 1 = ^ [a = - i i = 1 ] 



(58) 



(59) 



is equivalent to Weyl-squared model. However, it is more convenient to make use of 
the result of the previous section. Taking the identity satisfied by Weyl 2-form in four 
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dimensions, namely = 0, into account [19], the metric field equations (511 boil 
down to 

(*D*D-m 2 )R^ = *(R a A *C fia ) (60) 

where, despite the fact that curvature tensor in three and four dimensions have quite 
distinct properties, it acquires formal resemblance to the equations of its three di- 
mensional analogue, namely, the NMG theory equations in three dimensions, cf., 
Eqn. ( 1721 below. Note that, as has been found above, no QC terms involves in the 
trace of the equations in four dimensions and as a result, in the present case, one 
has R = 0. Consequently, ( f60b also entails Klein-Gordon type tensorial equation for 
Ricci 1 -form in four dimensions analogous to the NMG case studied in the following 
section. However note that it is not easy to see such a connection using equivalent 
lagrangian above. 

The linearization of the field equations of QC lagrangians supplemented with 
Einstein-Hilbert term in the metric perturbations reveals that in four dimensions gener- 
ically QC theories describe massless spin-2 modes, together with massive spin-2 
modes in addition to the massive spin-0 mode [?]. Although, QC terms renders 
four dimensional gravity renormalizable, unfortunately massive spin-2 excitations 
are ghost and there is a problem with unitarity. On the other hand, recently proposed 
Critical gravity model is obtained by the choice of constants a,b in as in the 
case (1) studied above together with addition of a particular cosmological term X * 1. 
The choice of + b = removes spin-0 mode whereas the cosmological constant 
with b + tV = effects the massive spin-2 modes that they become massless lfl2l . 

The linearization of field equations (f60l > around a vacuum solution in terms of 
metric perturbations immediately yield an equation with a fourth order linear dif- 
ferential operator which is of the form of a product of two second order differential 
operators. In addition, in three dimensions for NMG theory, the second order operator 
in (l60l can also be effectively factorized. The construction of Critical gravity in four 
dimensions also makes the models based on QC lagrangians studied the cases (2)-(3) 
more appealing for further study in n > 4 dimensions towards extension of critical 
gravity to higher dimensions [32|. 

7 Three dimensions 

Three dimensions can be considered as exceptional regarding the field equations for 
the gravitational theories that have second order trace. As is well-known, the Einstein- 
Hilbert action R * 1 does not have to a propagating degree of freedom and does not 
lead to a dynamical theory of gravity. On the other hand, there are two higher order 
theories which have desirable properties regarding the particle spectra and perturba- 
tive quantum theory based on these lagrangians. 

In this section, the classical aspects of Topologically massive and New massive 
gravity theories, which are third and fourth order in the metric components respec- 
tively, will be studied in the orthonormal coframe formulation provided above. It will 
explicitly be shown that the properties of the NMG field equations as a QC lagrangian 
are peculiar to three dimensions and in that the field equations NMG field equations 
are related to those of TMG theory in a special way under certain conditions. 
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The QC part of the lagrangian corresponds to the case (1) with Jz?j 2 [a = — \,b = 
1] and customarily denoted by J£k in three dimensions. In a more common form, 
NMG lagrangian explicitly reads 

J? NMG = aSe E H + = oQ a a A *e a V + —j{R a A *R a - 3 -R 2 * 1). (61) 

m L r m 8 

The sign of the Einstein-Hilbert relative to the QC terms has crucial effect on the par- 
ticle spectrum of the NMG lagrangian, thus the sign will be encoded in the constant 
(7. As before, the NMG field equations can be written down in the form 

-2m 2 O*G a + *T£+Dlg = (62) 

where the auxiliary 2-form X K P , that is to be used in *T^, takes the following equiv- 
alent forms 

xf = e a a (rP - -Re?) - eP a (R a - -R6 a ) 

8 8 
= £2 a P - ^R9 a P (63) 

in this case and using the general formulae (0 with TI^ = D*X^ leads to 

A# = 2*C" (64) 

where the result is a special case of the general expression (|29l ). The remarkable fact 
about d64l i is that all the fourth order terms are absorbed into the single expression 
involving the Hodge dual of Cotton 2-form. Since it involves the Hodge dual of the 
Cotton 2-form MNG equations are parity preserving unlike those of TMG In con- 
trast to the field equations of the QC lagrangian (l32l studied in the case (1) in n > 4 
dimensions, the fourth order lagrange multiplier term turns out to have the same form 
in three dimensions as well. In particular, recall that the simplified form of the la- 
grange multiplier term in n > 4 dimensions has been obtained by making use of the 
properties of the Weyl 2-form which in fact defined for dimensions n > 4. The result 
d64l > for the Lagrange multiplier forms for the NMG theory [20 1 also implies that the 
NMG field equations can be related to TMG equations both of them involve Cotton 2- 
form. Moreover, the Schouten 1-form can be regarded as a kind of potential because 
of the relation C a = DL a . 

Note that in both TMG and NMG field equations, higher order terms involve 
Cotton 2-form, or equivalently co variant exterior derivatives of the Schouten 1-form. 
Moreover, as a consequence of the trace-free property of the Cotton tensor, cosmolog- 
ically extended TMG field equations has the property that the trace of the metric field 
equations are second order as in the case of NMG. As will explicitly be shown below, 
such similarities implies closer connection between the mathematical structures of 
the two theories. 

In three dimensions, it is well-known that there is no Weyl 2-form and the con- 
formal properties of three dimensional metrics are encoded in Cotton 2-form, see, 
for example, ll23l . Consequently, it is evident from the NMG equations (l62l that for 
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conformally flat solutions, the order of the NMG equations are in fact second as a 
result of the fact that — for conformally flat solutions. 

The NMG field equations written in the from (l62l l further implies that it is possible 
to relate them to TMG field equations when the field equations for two theories are 
written in terms of Schouten 1-form. In order to provide such a connection between 
NMG and TMG, it is first convenient to recall that TMG field equations follow from 
gravitational Chern-Simons action supplemented with Einstein-Hilbert action iflOl 

EE] 

1 2 

■&TMG = °R*l + 2-(fl»p At/ffl^ + ^0) a p A A C0^ a ) (65) 

where a multiplicative constant cr is introduced for convenience. In terms of differ- 
ential forms, the vacuum TMG field equation that follows from this lagrangian reads 

a*G a + ^C a =0 (66) 

where fl ^ = ±QP a and X a = ~^L a for the Chern-Simons part in (|65). The trace 
of the vacuum TMG field equations imply that R = as a result of the fact that the 
Cotton tensor has vanishing trace i a C a = 0. Now, it will be convenient to rewrite (16*6l l 
in terms of Schouten 1-form as 

{*D-oii)L a = (67) 

by taking the Hodge dual of 1661 . Then by applying the operator *D to d67b and using 
d67l > in the resulting expression, one finds 

{*D*D-o 2 ll 2 )L a = (68) 

which are fourth order in the metric components relative to a coordinate basis. Note 
also that for R = it is possible to replace L a with R a in 068) . By linearizing the 
full nonlinear Eqn. d68l by the replacement of L a with the corresponding linearized 
Ricci tensor and at the same time by the replacement of the sequence of the covariant 
exterior derivatives and Hodge duals *D* D with three dimensional Laplacian of 
the Minkowski background, one finds that the linearized curvature (in the traceless 
transverse gauge) represents an excitation spin-2 particle with mass jj. [ 10] . 

As will explicitly be shown below, the form of the field equations, in fact, suggests 
that NMG can also be cast into a form analogous to (l68l . In passing note that inclusion 
of a cosmological constant results in a constant scalar curvature and consequently a 
covariantly constant source term to (|68V The cosmologically extended TMG equation 
then becomes 

{*D*D-o 2 ll 2 )L a =q a (69) 

where q a = j(jA\p. 2 9 a is covariantly constant vector-valued 1-form: Dq a = 0. In the 
form given above, the highest derivative term in the squared TMG equations contain 
the term D*C a which is also present in the NMG equations d62] l. 
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Now, one can work from NMG side towards (l6?t or (169V By using the definition 
C a = DL a and also expressing the Einstein forms in terms of Schouten 1 -forms and 
scalar curvature, NMG field equations (l62t can be rewritten in the form 

{*D*D + om 2 )L a = ^<jm 2 R9 a + ^Tg (70) 

The source term on the right hand side can further be simplified by using the trace 
of the field equations. Explicitly, by wedging d62l from the left by the basis coframe 
form 9jj one finds 

2m 2 (jR*l+£f K = 0. (71) 
Consequently, the NMG field equations take the form 

{*D*D + om 2 )L a = l ~(Tg- l -T K e a ) (72) 

where Tk = is the trace. It is worth emphasizing that, up to this point, neither any 
simplifying assumption nor any approximation has been introduced into the deriva- 
tion of NMG equation d72b and it holds in full generality. The formal similarity be- 
tween (172b and d68l i becomes even more refined in the cases where 

(i) The source 1 -form in (l72t vanishes identically, 

(ii) It is proportional to Schouten 1-form, 

\{Tg- l -T K e a )=kL a . (73) 

In the latter case, the NMG field equations take the form 

(*D*D-m 2 )L a =0. (74) 

where the shifted mass has then been defined by — m' 2 = am 2 —X. Thus, as a re- 
sult, if a solution of NMG satisfies the algebraic constraint ( l73l . also satisfies the 
Klein-Gordon type Equation ( f74b and therefore it satisfies the "square" of the TMG 
equations provided that the mass parameters of the theories are simply related by 
jll = m'. 

The foregoing analysis can easily be extended to the case of TMG and/or NMG 
extended by respective cosmological constants A\ and A%. For instance, for the latter, 
the NMG field equations d72b extended by cosmological constant (via the additional 
term — 2A 2 * 1 in the NMG action) take the form 

{*D*D + om 2 )L a = ±(T£ - l -T K e a +A 2 m 2 e a ). (75) 

An immediate consequence of Eqns. (l74l > and/or (l69li-(l75ll is that with the alge- 
braic constraints specified above satisfied, the TMG and NMG theories admit com- 
mon set of solutions. Consequently, with the algebraic relation between the mass pa- 
rameters of the two theories and the condition relating expressed in terms Schouten 
1-form and QC energy momentum 1-form T%, TMG field equations are in fact can 
be considered as the exact "square-root" of the NMG field equations. This corre- 
spondence has previously been introduced in [36 1 in a more general context and used 
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to find the solutions of NMG field equations using the field equations of TMG in 
a series of papers ll36l . The presentation given here in terms of exterior calculus of 
forms allows one to identify the Dirac type operator introduced in 1 36 1 simply as 
= *D in terms of geometrical operators, namely, the covariant exterior derivative 
in composition with the Hodge dual relative to an orthonormal coframe. 

Now it is a convenient point to comment briefly on the equations for General 
massive gravity (GMG) in the framework above. The GMG theory follows from the 
lagrangian which is a combination of those of TMG and NMG theories and thus 
has two distinct mass parameters [11|. TMG and NMG theories can be recovered 
by taking appropriate limits of the two mass parameters of GMG. It is convenient to 
consider directly the field equations for GMG which are of the form 

a*G a + -C a + \(D*C a + *Tg)=0. (76) 
By introducing the mass parameters for GMG by redefinitions m 2 = m+m_ and 

2 

^- = m+ it is easy to show that the GMG equations take the form 

[*D * D - (m+ +m-)*D + m + mJ\L a = (77) 

in terms of the new mass parameters with the QC source term Tjj- is defined as the 
same as that of the NMG theory. Note that (1771 1 is not a linearized equation in metric 
perturbation and holds in full generality. Thus, it is not possible to write the differen- 
tial operator on the left hand side in a factorized form. On the other hand, the equation 
provides a convenient expression with regard to linear approximation in the metric 
perturbation around a flat background. Explicitly, the leading term on the left hand 
side of d771 > is of the form 

[*D *D (m + + m_) *D + m + mJ\L a = [*d *d— (m + + m) *d + m + mJ\L a + ... 

« (*d - m+)(*d - m-)Lf in + ... (78) 

where the Hodge duals on the right hand side on the second line are assumed to be 
approximated by the Hodge duals of the flat background. The factorized differential 
operators of the flat background act on linearized Schouten 1 -form Lf. . As the GMG 
calculations elegantly illustrate, the mathematical formalism presented allows to ob- 
tain the linear terms in a factorized form and in fact these considerations apply to the 
other QC field equations above. 

Finally, recall that a constant curvature spacetime, with curvature 2-form, 

^ = -—^—-0^, (79) 

(« -!)(«- 2) 

satisfies the Einstein field equations *G a = A * 9 a . For the general lagrangian dm 
supplemented with Einstein-Hilbert action and a cosmological constant Aq, i.e., for 
the lagrangian of the form 



W = - (n aP a * e a P - 2A Q * l ) + _s?/ 



(n) 



(80) 
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the Einstein metric of (l79l l produces 

l2 4(«-4) 



*T, L 



(11-2)2 



(II -1) 



*0 a 



(81) 



which can be regarded as a kind of cosmological term. In addition, one has DX, a — 
identically. Eventually, with the maximally symmetric metrics, the vacuum field 
equations for ( f80b reduce to 



Aq-A- 



2(n-4) 
K-(n-2) 2 



(b-1) 



Thus, if the cosmological constant satisfies the equation 

an + b 



(n-1) 



*e a = o. 



o 



(82) 



(83) 



then, d79i l also solves the vacuum equations of the QC model ( f80b [35|. However, 
now there are two distinct vacuum solutions corresponding to the roots of Eqn. (1831 1. 
Thus, generic QC field equations can be linearized around the maximally symmetric 
vacua as well as flat vacuum background. 



8 Concluding remarks 

In the above study of QC gravity equations, a concise form of the general QC gravity 
equations < fT~5T > relative to orhonormal coframe are derived. These provide a different 
point of view and further insight into the QC field equations. They can be cast into 
various other forms that may be useful for some other purpose s as well. In this regard, 
note that d49b can be rewritten in the form 

£l a P = j^wf (84) 

k=\ 

which shows explicitly the irreducible decomposition of the curvature 2-form [34|. 
W"^ is traceless fourth rank part W"^ = C a $ . W"^ is traceless Ricci part W"^ = 
- r - L ^(9 a ASl 3 - el 3 AS a ) where S a = R a - ±R9 a is the traceless Ricci 1-form and 

finally, W" P = -^Y)RO afi is the trace part. Then, using this decomposition of the 
curvature 2-form, the second equality in (l50t can be written in the form 

3 

Q a P A *i2 a j3 = 52 W k P ^*W kap . (85) 

k=l 

Moreover, the irreducible parts are orthogonal in the sense that n a P A*W ka B = 
Wk a $ A *Wk a p and Wk a P A *Wj a p x Sjk- Now, it is easy to formulate the QC models 
by employing the lagrangian densities of the form Wk a ^ A *Wj a p. The correspond- 
ing bivector valued 2-form then takes the form X"^ — 2Wk a ^- Therefore, in terms of 
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these, the field equations ( TT3T > directly expresses the field equations in terms of the 
contractions of the irreducible parts of the curvature. Consequently, (fT5T l can be put 
in a convenient to form to write the general QC field equations relative to the null 
coframe of Newman-Penrose in terms of spinor components 1137 113 81 . 

Apparently, the construction of QC lagrangians having second order trace equa- 
tions in terms of the bilinear combinations considered above requires the coupling 
constants (free parameters) that depend crucially on the number of spacetime dimen- 
sions. As a consequence, a QC lagrangian density having this property for a given 
spacetime dimensions does not generalize to other dimensions. Regarding the trace 
property, each spacetime dimension has to be considered separately. In four dimen- 
sions, for example, the most general QC gravitational lagrangian that has second or- 
der trace is of the form (l3Zt with particular values of the coupling constants a = — \b. 
However, for n > 4, the most general QC lagrangian contains all the three terms con- 
sidered in the cases (l)-(3) above that cannot be eliminated by adding a topological 
term as in four dimensions. Therefore, there are three independent linear combina- 
tions of QC lagrangians each of which has the desired property of second order trace 
and consequently, the most general QC lagrangian contains three free parameters. In 
dimensions n > 4, the QC lagrangians derived in cases (2) and (3) above deserve fur- 
ther scrutiny with regard to, for example, their particle spectra and classical solutions 
as well. 

One of the original results of the paper is that, the orthonormal coframe formu- 
lation of the field equations allows the lagrangians of the form ( f58T l to be cast into 
Klein-Gordon type tensorial equation of the general form (l72l in any dimensions 
n > 3. Three dimensions seem to have peculiar properties that has no analogues in 
higher dimensions. TMG and NMG field equations are related in a special way, in 
that, with certain conditions satisfied, the former, which can be written as a Klein- 
Gordon type equation considered as the square-root of the latter. The direct approach 
of solving QC gravity equations starting with a suitable metrical ansatze is difficult, 
however, the study of coframe formulation of QC field equations may provide insight 
into the QC field equations. In this regard, further attention is required to provide con- 
crete examples of class of metrics that satisfy the conditions either (i) or (ii) which 
makes the TMG equations square root of the NMG equations. 

Another original result is the construction of general QC models having second 
order trace in n > 4 dimensions with three free parameters. Consequently, as can 
be observed from Eqn. (146V some these QC models involve contraction of Riemann 
tensor with itself and they have no analogues in n < 4 dimensions. For these models, 
the field equations cannot be written in terms of Schouten tensor as opposed to the 
well-known QC models involving linear combinations of Ricci-squared and scalar 
curvature-squared terms in n > 3 dimensions. 

In four dimensions, it may provide further insight for the field equations of the 
Weyl-squared model if the field Eqn. (Io*0l i can be rewritten in terms of the Lanczos 
potential for the Weyl 2-form [39|. This, however in the present framework can nat- 
urally be achieved by making use of the Lanczos potential formulated relative to an 
orthonormal coframe and this requires separate investigation. The field equations ex- 
pressed in terms of the exterior differential forms provide further insight, at the least, 
into the NMG field equations in relation to TMG theory. Written down explicitly in 
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terms of differential forms, and also in terms of Schouten 1-form, the novel descrip- 
tion of the NMG field equations ll36l amounts to taking Hodge dual of the metric 
field equations and subsequently regrouping the QC terms conveniently. This result 
follows from the fact that relative to orthonormal coframe all the fourth order terms in 
the field equations can concisely be expressed in terms of Cotton 2-form for the NMG 
theory. In contrast, the novel description of NMG in terms of traceless Ricci tensor 
given in [ 36 1 constructed by using identities satisfied by three dimensional curvature 
tensor. 

Finally, it is interesting to note that the second order covariant differential operator 
D * D* which comes into play by construction of the NMG field equations relative 
to an orthonormal coframe also shows up in a quite different context of introducing 
mass to graviton in other gravitational models. In [40 1, using the first order formalism 
ll34l in three dimensions and in the context topological Riemann-Cartan type theory 
generalizing TMG, the basis coframe 1 -forms are shown to satisfy a Klein-Gordon 
type equation of the same form as ( f74l i in terms of the differential operator *D * D =p 
D*D*. 
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